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ABSTRACT 



Context. Magnetic fields at the surface of a few early-type stars have been directly detected. These fields have magnitudes between a 
few hundred G up to a few kG. In one case, evidence of magnetic braking has been found. 

Aims. We investigate the effects of magnetic braking on the evolution of rotating (u m i=200 km s _1 ) 10 M Q stellar models at solar 
metallicity during the main-sequence (MS) phase. 

Methods. The magnetic braking process is included in our stellar models according to the formalism deduced from 2D MHD simu- 
lations of magnetic wind confinement by ud-Doula and co-workers. Various assumptions are made regarding both the magnitude of 
the magnetic field and of the efficiency of the angular momentum transport mechanisms in the stellar interior. 

Results. When magnetic braking occurs in models with differential rotation, a strong and rapid mixing is obtained at the surface 
accompanied by a rapid decrease in the surface velocity. Such a process might account for some MS stars showing strong mixing and 
low surface velocities. When solid-body rotation is imposed in the interior, the star is slowed down so rapidly that surface enrichments 
are smaller than in similar models with no magnetic braking. In both kinds of models (differentially or uniformly rotating), magnetic 
braking due to a field of a few 100 G significantly reduces the angular momentum of the core during the MS phase. This reduction is 
much greater in solid-body rotating models. 
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1. Introduction 

Many indirect observations indica te that magnetic field s should 
be pr esent around massive stars ( He nrichs et alj l200a : iHubrigi 

2008) , and recently a few direct detection s have been made via 
the Z eeman effect (see the recent review bv lDonati & Land street 

2009) . For instance, direct measurements have been obtained for 
a handful of O stars, among t hemfl 1 Ori CandHD 19 1612, with 
values of a few hundred G (Donatieta 02TM 12006a). This last 
star is a member of the class of Of?p stars, of which only 5 are 
known in the Galaxy. Of these, three have bee n observed in deta il 
and all three have detectable magnetic fields (IWade et alJl2010h . 

Among early B-type stars, a m agnetic field has be en dis- 
', covered in the B0.2V star r Sco dDonati et all l2006bl ~ 0.5 
kG) and in three /3 Ceph eid stars: £' CMA (iHubrig et al.1l2006l 
B0.7IV, 300 G), p Cep dHenrichs et alj I200QL B2III. less than 
100 G) and V2052 Qph dNemer et al.N2003b[ B2IV , less than 
100 G). INeiner etail d2003al) and lHubrig et all d2006l) have dis- 
covered magnetic fields of the order of a few hundred G in 
slowly pulsating B-type stars. Magnetic fields have been de- 
tected around a few Be star s dYudin et al.l l2009t IHubrig et al.l 
l2007t INeiner & Hubertl 120051) . and these fields are of the order 
of 100 G or less. 

The origin of these magnetic fields is still unknown. It 
might be fossil fields (e.g. the spectral characteristics of Of?p 
stars are indicative of organi zed magnetic field s, most likely 
of a fossil origin according to lWade et al.l l2010). or fields pro- 
duced through a dy namo mechanism. Recent simulations by 
ICantiello et al]d2010b of subsurface convective zone in massive 
stars show dynamo-generated magnetic fields of the order of 
one kG. According to these authors, these magnetic fields might 
reach the surface of OB stars. 



In the Sun, magnetic braking results from solar wind material 
following the magnetic field lines that extend well beyond the 
stellar surface. This coupling exerts a torque on the surface lay- 
ers of the Sun, and this slows down its rotation. Could such a pro- 
cess also be active in m assive stars showing a sufficiently strong 
surface magnetic field? iTownsend et alJ d2010h has recently dis- 
covered that cr Ori E (B2Vpe, ~ 10 kG) is undergoing rotational 
braking. The spin-down time of 1.34 Myr is in good agreement 
with theoretical predictions based on magnetohydrodynamical 
simulations of angular momentum loss from a magnetized line- 
driven wind. This gives some support to the hypothesis that at 
least a few massive stars may indeed suffer magnetic braking. 

Such an effect has for the moment never been included in 
massive star models; however, by modifying the internal distri- 
bution of Q, the angular velocity, magnetic braking can signifi- 
cantly change the mixing of the elements inside the star, as well 
as the evolution of its angular momentum content. 

Although rotating models have improved the agreement be- 
tween models and theory, some points remain to be clarified. A 
small subset of stars exhibit surface properties, such as low sur- 
face velocities and strong surface enrichments, which, when the 
star is not a giant or a supergiant, cannot be explained by cur- 
rent rotati ng stel lar models fo r sing le stars (see discussions in 
lBrottetalj|2009T; iHunter et alj|2009l) . Although the fraction of 
these s tars is small and s ome may be stars at the end of the MS 
phase dMaeder et al.ll2009h . it is worthwhile investigating what 
the physical cause of this behaviour could be. 

Theoretical models with rotation also predict rotation veloc- 
iti es that are too high for th e young pulsars (see the discussion 
in iHeger et al J 12004 l2005h . Thus it does appear interesting to 
study any effects that can remove angular momentum from the 
star. Can magnetic braking be an interesting explanation for the 
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Fig. 1. Evolution of the N/C ratio (in number) as a function of 
the gravity at the surface of 10 M stellar models with Vj n ; = 200 
km s _1 . The value of the magnetic field is indicated. Continuous 
(dashed) lines are for models with interior differential (solid 
body) rotation. In the shaded areas, models with differential ro- 
tation predict surface velocities inferior to 50 km s _I during the 
MS phase. The full circle and square correspond to HD 16582 
and HD 3360 (£ Cas) respectively (see text for references). 

strongly mixed, slow rotators? Can it help in reducing the angu- 
lar momentum of the core? In this first paper, we want to study 
these questions by presenting a first series of computations ac- 
counting for this effect. In section 2 we present the formalism 
we used to account for the magnetic braking effect, as well as 
the physical ingredients of the models. The results are discussed 
in sections 3 and 4 and conclusions are given in section 5. 

2. Physics of the stellar models 

The formalism of the magnetic braking law used here follows 
theoretical developme nts first made for the Sun, starting with 
Web er & Davisl (1 19671) who used an idealized monopole field to 
model the angular momentum loss in the solar wind J. They 
found that j = 2/3MClR 2 A , where M is the mass loss rate and 
Ra the Alfven radius, defined as the point where the ratio be- 
tween the magnetic field energy d ensity and the kinetic energy 
densi ty of the wind is e qual to 1 . lud-Doula & Owoc ki (2002) 
and ud-Doula et al ] (120081) have examined the angular momen- 
tum loss from magnetic hot stars with a line-driven stellar wind 
and a rotation-aligned dipole magnetic field using 2-D numerical 
MHD simulations. They find that the total angular momentum 
loss follows the expression given above, but with a smaller 
than in the monopole case. They find that 

— = -MQRI[Q.29 + (77, + 0.25) 1/4 ] 2 , (1) 
df 3 

where ri„ is the m agnetic confinement parameter 
(lud-Doula & Owockil |2002|) defined by 77, = B^Rl/Mv^, 
where R„ is the stellar radius, B eq , the magnetic field at the 
equator and at the surface of the star, and Voo the wind velocity 
at infinity. Very interestingly, when applied to the case of cr 
Ori E, Eq.Q] gives a spin-down timescale of the order of 1 Myr 
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Fig. 2. Same as Fig. [T]for the evolution of the abundance of B at 
the surface. 



dud -Doula et alJ l2009h . quite in agreement with the o bserved 
spin-down timescale for this star dTownsend et al.ll2010h . 

We have imp l ement ed in the Geneva evolution code 
(Eggenber ger et all 120081) Eq. Q] for the loss of the angular 
momentum. Various values of B eq are used between 0.1 and 
10 kG. In the present exploratory work, we suppose that the 
magnetic field keeps a constant value during the MS phase. 
We computed the evolution of 10 M stellar models at so- 
lar metallicity with v\„j = 200 km s _1 . In the context of shel- 
lular rotation dZahn| 1992), the vertical transport of chemicals 
through the combined action of vertical advection and strong 
horiz ontal diffusion can be d escribed as a pure diffusive pro- 
cess dChabover & Zahnlll992l) . The advective transport is then 
replaced by a diffusive term, with an effective diffusion coef- 
ficient D e ff = |ri/(r)| 2 /30Z)h, where r is the characteristic ra- 
dius of the isobar, U(r) the vertical componen t of the merid- 
ional circulation velocity dMaeder & Zahnll 998) and the dif- 
fusion coefficient associated to horizontal turbulence as given by 
IZahnl dl992l) . The vertical transport of chemical elements then 
obeys a diffusion equation, which in addition to this macroscopic 
transport, also accounts for (vertical) turbulent transport with the 
same coefficient Z) s h ear as for the transport of angular momentum 
dMaeder & MevnefcOOlh . 

The effects of magnetic braking depend a lot on how angular 
momentum, hence the chemical species, is transported inside the 
star. Here we consider two limiting cases: 

- Differential rotation: the angular momentum transport is 
driven by the meridional currents and the shear instabilities. 
These two effects are not efficient enough to produce solid 
body rotation. A moderate differential rotation is therefore 
present during the whole MS phase, and the mixing of the 
chemical species is mainly due to shear instabilities. 

- Solid body rotation: when the transport of the angular mo- 
mentum is very efficient, then solid body rotation is main- 
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Fig. 3. Same as Fig. Q]for the evolution of the surface velocity at Fig. 4. Variation in the specific angular momentum as a function 
the equator. of the Lagrangian mass inside the cores of different models. 



tained during the whole MS phaseQ. The chemical species 
are transported by the meridional currents. 

At the present time it is not known which is the most frequent 
situation in nature, so we explore these two possibilities, which 
represent some limiting cases hopefully framing the real world. 



3. Models with magnetic braking and interior 
differential rotation 

In Figs. Q] [2] and[3](look at the continuous lines) are shown the 
evolutions of the abundances and of the (equatorial) velocities at 
the surface of models computed with different values of the mag- 
netic field. We see that, compared with the model without mag- 
netic braking, models with magnetic braking present stronger 
changes of the surface abundances. This comes from the fact 
that magnetic braking creates strong differential rotation in the 
outer layers and thus triggers mixing by shear. At the same time 
the surface velocity rapidly decreases. 

More quantitatively, we see that, at 100 G, the effects on 
N/C remain modest (increase by about 20% at the end of the 
MS phase with respect to the model without magnetic braking) 
but are already significant for the surface velocity, which is de- 
creased from 168 to 115 km s _I at the blue turnoff of the MS 
band. For a 500 G magnetic field, the effects are very strong: in- 
crease by more than a factor 2.2 of the N/C ratio at the end of the 
MS phase. The surface velocity reaches values around 19 km s _1 
at the end of the MS phase. Surface abundances of boron provide 
a very sensitive test of how mixing occurs in the outer layers of 
the star. Indeed, that element already begins to be destroyed at 
temperatures between 5-6 10 6 K. From Fig. [2] we see that in 
models with magnetic braking, boron is much more rapidly de- 
pleted at the surface. 

Stars in the hatched zone of Figs.Q]and[2]present strong signs 
of mixing during the main sequence and have surface veloci- 
ties inferior to 50 km s . Typically, one would expect that at 

1 We do not invoke any specific mechanism here. Solid-body rotation 
might be due to a fossil magnetic field fo r instance or to a dynamo 
driven by shear such as the one proposed by Spruit ( 2 0021) . 



least some LMC stars in group 2 of Fig. 5 in lHunter et al.l |2009) 
would be in this regi o n. We defer a detailed comparison with the 
data of iHunter et all (120091) to a later paper where models with 
the LMC metallicity will be computed. In Figs.Q~|[2] and [3] we 
have indicated the positions of the two known stars having initial 
masses between 9 and 1 1 M and showing boron depletion by 
m ore than an order of mag nitude (data taken from Tables 3 and 4 
of iFrischknecht et al]|2010i see detailed references therein). We 
see that current models with magnetic braking although not pro- 
viding a perfectly consistent fit to the dat£0, appear to be promis- 
ing for explaining stars with strong signs of mixing, low surface 
velocities, and high gravities. It is interesting to mention that the 
star corresponding to the square (<T Cas) has a detect ed magnetic 
nolar field of 335f+120-65 N ) G (iNeiner et al.ll2003ah . 

Can this process extract angular momentum from the core 
region? In Fig. [4] the specific angular momentum inside the cen- 
tral 2 M Q at the end of the MS phase is shown for various models. 
We recall that the reduction between the angular momentum of 
the core on the ZAMS and at the presupernova stage occurs for a 
significant part during the MS phase. Typically about 70% of the 
total angular momentum lost by the core during the whole stellar 
lifetim e is lost during the MS phase (see Fig. 7 in lHirschi et al.l 
2004). In models with differential rotation, a magnetic field of 
about 100 G reduces the specific angular momentum by about 
25% (0.1 dex) with respect to the values obtained in the model 
without magnetic braking. When a magnetic field between 500 
G and 1 kG is considered, j is on average decreased by slightly 
more that 0.3 dex, which means by a factor between 2 and 2.2 
with respect to the model with no magnetic braking. Thus mag- 
netic braking in differentially rotating models has an impact on 
the evolution of core angular momentum content. 
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4. Models with magnetic braking and interior solid 
body rotation 

The case of solid body rotation is quite different from the one of 
differential rotation as can be seen in Figs. Q] [2] [3] and[4] This is 
expected since the physical processes responsible for the trans- 
port mechanisms are quite different. For the angular momentum, 
a strong coupling is imposed in the model. The chemical ele- 
ments are no longer transported by shear turbulence as in the 
previous case, but by meridional currents. 

Whenever meridional currents are responsible for the trans- 
port of the chemical species, the key factor governing chemical 
element mixing is Q, which decreases rapidly in the whole in- 
terior when magnetic braking is exerted at the surface. When 
magnetic braking is accounted for, mixing is much weaker than 
in the case without magnetic braking. For a magnetic field equal 
to or larger than 0.5 kG, no nitrogen enrichment is predicted by 
the models during the whole MS phase. Boron depletions are not 
as strong in models with magnetic braking. 

From Fig.[3j we can see that for models with magnetic brak- 
ing and solid body rotation, the surface rotation decreases much 
more rapidly than in the case of differential rotation (for a given 
value of B eq ). Thus such models would predict stars at the end 
of the MS phase with low surface velocity and weak or even 
nonexistent surface enrichments. Such stars can also result from 
the evolution of progenitors starting with low initial rotational 
velocities. However, the internal structure at the end of the MS 
will be different if the star has ever been a slow rotator or has 
come from a slowed down, initially rapid rotators. Typically the 
star with a fast-rotating progenitor will have a larger core. 

Solid body rotation alone (without magnetic braking) is al- 
ready an interesting way to evacuate the angular mo mentum en- 
closed in the core. For instance, it has been shown by Hege r et al.l 
(20051 120041) . who computed models where solid body rotation 
is imposed duri ng the MS p hase by the action of the Tayler- 
Spruit dynamo dSpruill 12002). that these models could produce 
neutron stars with rotational velocities at birth in the upper range 
values given by observations. As can be seen from Fig.|4](lower 
panel), magnetic braking with solid body rotation has a dramatic 
effect on the angular momentum content of the core. We see that, 
already with a value of B ecj = 100 G, j is decreased by more than 
one dex with respect to the corresponding solid-body rotation 
with no magnetic braking. For magnetic fields of 1 kG, j is de- 
creased by more than 4 orders of magnitude! In that last model, 
the specific angular momentum in the core at the end of the MS 
phase has the same order of magnitude as the one required for 
explaining the lon g-period (150 ms ) young pulsar PSR B1509- 
58 (see Table 3 in lHeger et aLEool . 

5. Conclusions 

The main results in this paper are the following. Models with 
differential rotation and magnetic braking may produce strongly 
mixed stars with low surface v elocities during the MS phase. 
Some of the stars discussed by Hunt er et al.l (I2009L see their 
Group 2 stars in their Fig. 5) as well as some sta rs showing 
very strong boron depletion with no N-enhancement dVenn et al.l 
l2002t iMendel et al.ll2006t iFrischknecht et al.ll2010h . might be 
stars that suffered magnetic braking. Models with solid-body ro- 
tation and magnetic braking produce stars that at the end of the 
MS phase have low surface velocities and none or very weak 

2 All our models start their evolution with a velocity of 200 km s -1 , 
which may not be adequate for these stars. 



changes in the surface abundances. Magnetic braking reduces 
the angular momentum content of the core, even for moderate 
values of the magnetic field when solid body rotation is consid- 
ered. 

Present models are of course quite preliminary and have to 
be improved in the future. For instance, the coexistence of a 
large-scale/dipolar magnetic field (as adopted here) with inter- 
nal differential rotation may be problematic. If the field is of fos- 
sil origin, it will likely pervade the whole star and may induce 
strong coupling, thus reducing or even removing any differen- 
tial rotation. In that case, our models with solid-body rotation 
would apply. In these models, magnetic braking is so efficient 
that the slow down occurs before any significant changes in the 
surface abundances. If the surface magnetic fiel ds are produced 
by a d ynamo (see for instance the suggestion bv lCantiello et al.l 
2010), differential rotation may still be present in the interior, 
and changes in the surface abundances are quite rapid. However, 
one may question whether it is justified to assume a constant 
value for B eq as we did here. More quantitative estimates need 
detailed understanding of the dynamo process and of its de- 
pendence on rotation. Finally let us stress that in the Hunter 
diagram most stars follow the predictions of rotating models 
(Mae der et al1 [2009) without magnetic braking. This is consis- 
tent with the weak magnetic fields of OB stars. 
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